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iBSTR.AGT 


The prohlem of Pil«l symbol deteotion over fibre optical 
digital data transmission system is considered. A sub- 
optimum solution is obtained under the assumptions that 
all past symbols have been decoded correctly j the electron 
emission times are quantized and signal to noise ratio is 
high. An upper bound on probability of error is numerically 
calculated using Gram-Oharlier series expansion. The per- 
form.a.nce is compared with that of an 'Integrate and Dump’ 
receiver, 'One Baud' receiver (when future symbol ISI is 
neglected) and an 'Ideal' receiver (in which future symbols 
are assumed to be known) . Results are also compared with 
simulation studies. Timing error effects and effect of 
sampling rate are also considered. 
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CHAPTER 1 


INTRODUCTION 

1 large part of traditional communication system 
practice is directed to detecting and processing electrical 
signals received over diverse time-continuous communication 
channels. Tue physical realization of each of these tra- 
ditional systems has led to mathematical treatment designed 
to handle problems such as linear distortion or fading, 
which were peculiar to one or even several such systems. 

But the principal concern of mathematical treatment of 
these time-continuous channels has been the ubiq.uitous 
additive Gaussian noise. In fact, it would be fair to say 
that much of the structure of the mathematical treatment 
used has been dictated by the mathematical properties of 
this noise. In the absence of noise, many problems would 
immediately degenerate, at least theoretically, to 
situations of perfect detection, infinite capacity etc. 

The consideration of some promising optical communi- 
cation systems seems to alter the above picture. Here, we 
have in mind the transmission of information by way of 
light pulses propagating through an optical fibre and 
subsequently detected by a photo-detector that converts 
electromagnetic energy in the fibre to electrical signals 
in a circuit. We immediately note certain features which 
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this problem has in common with the traditional problems. 
For one thing, the fibre can delay, attenuate or spread' 
the transmitted pulses. For another, the electrical signal 
after photo-detection may be corrupted by additive Gaussian 
noise. Yet there is another fundamental impairment. The 
received intensity of the electromagnetic signal that 
propagates in the fibre (which acts as a waveguide) is, 
under practical conditions, sufficiently weak, ind any 
effective scheme of detection of the signal must take 
into accoimt the quantum nature of the electromagnetic 
disturbance. In other words, the detection mechanism must 
be based upon photon counting. Here, a new element enters 
the problem: phot on counting is subject to statistical 
fluctuations. In the quantum case, a signal un corrupted 
by any external disturbance still carries with it, its own 
'noise', as it were, which is not Gaussian. This noise, 
described as an electron counting process, is a time- 
varying Poisson process whose intensity (or rate) function 
pN(t) varies in direct proportion to the inf ormation-bearing 
pulse train. 

This thesis is devoted to a study of determining a 
receiver structure for the above mentioned situs. tion and 
evaluating its performance. 

1.1 ORGMISATIOH CF THE THESIS 

The thesis comprises 5 chapters. Chapter P^after 



introducing a mathematical model for the fibre optical 
communication system, briefly describes the various tech- 
niques proposed in the relevant literature, for data 
transmission over the optical fibre. 

In Chapter 3 a receiver structure is derived and a 
hardware block schematic for the receiver is presented. 

This is followed by a discussion on the complexity of 
the proposed receiver. 

In Chapter 4 the performance of the proposed equalizer 
is analysed and it is compared with that of an 'Integrate 
and Dump' receiver, 'One Baud' receiver (where future 
symbol ISI is neglected) and 'Ideal 'receiver (in which 
future symbols are also assumed to be known) . Some simula- 
tion results are also presented. 

Chapter 3 concludes this thesis, describing some short- 
comings of the present work, and suggests approaches for 
further work to refine the receiver described herein. 



CHiiPTBR 2 


A MODEL FOR DIGITAL FIBRE OFTIC BTSTBM AND 
LITERATURE REVIEW 


In this chapterjcf irst we present a mathematical model 
of a typical fibre optical communication system. Then we 
formulate the optical detection problem. After this we 
describe the previous work done to solve this problem. 

2.1 A MODEL FOR DIGITAL FIBRE OPTICAL COMBTONI CATION SYSTEM 

Here we consider the binary pulse amplitude modulated 
(PAM) transmission system. Here the choice between a 'one' 
or a 'zero' is translated into the presence or absence of 
a short burst of optical power (light) in the fibre. To 
understand this in more detail, we shall trace the passage 
of a single pulse through the mathematical model of our 
system (see Fig. 2.1). In case of an 'one' being trans- 
mitted, an electric pulse (a pulse of duration T) turns 
on a 'flash light', which in this case is a laser or light- 
emitting diode, and electromagnetic energy is sent into the 
transmission medium (optical fibre). In Fig. 2.1 we have 
shown how this pulse spreads as it travels along the 
fibre. Fig. 2.2 shows some of the typical pulse responses 
of the fibre optical system. These curves are taken from 
Cartledge [3]. There are two courses open to us when this 
pulse reaches the other end of fibre. First, we can process 
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the received electromagnetic field dictated by quantum 
mechanical considerations [4]. But such schemes of 
processing are complex. Instead, we adopt a practical 
approach and assume that a direct detection is used to 
convert optical energy to an electrical signal. This is 
accomplished by using a photodetector prior to any signal 
processing. The electrons will be emitted due to the 
electromagnetic energy falling on the detector. Exactly 
when in time the electrons are emitted is random and is 
the Poisson process mentioned earlier. The probability of 
receiving a count between time t and t+dt is proportional 
to p(t)dt, whe^e, owing to effects in the fibre, p(t) is a 
distorted and attenuated version of the transmitted pulse 
(see Fig. 2.2). In practice a background of counts also | 
exists. This is called the dark current and is modelled by : 

introducing a constant additive intensity function before 

* !. 

the detector, although some of these counts can originate ^ 

! 

in the physical detector itself. Typically, transmitted J 
power and transmission loss (fibre length) are adjusted so | 
that, of the order of one to two hundred photons per pulse j 

f 

are, on the average, detected. The dark current contribute^ 
about 1 to 5 percent of counts. 

To transmit a zero, we simply do not turn on the 
transmitting power, and detector only registers counts [ 

resulting from the dark current. '■ | 
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We have been loosely speaking of the output of the 
photo-detector as 'counts’* The actual current at the 
output of this device caused by a photon is a wide-band 
pulse g.w(t) (very narrow compared to T, a delta function 
in the limit), where g is an integer valued random varia- 
ble or g = 1 , depending on whether or not an avalanche 
diode issued. The electric current at the output of the 
photo-detector is further distorted by Gaussian noise 
whose effects are of lesser importance when an avalanche 
diode is used but not for the case when g = 1 . The pulse 
energy-to-noise ratio, is typically -20 dB, 

In other words, the response of the photodetector to an 
individual photon is masked by the additive back ground 
noise. This is in contrast to the photomultiplier and 

avalanche photodiode whose mean gain g and mean square 
-2 

gain g can be rather large. For these devices the average 
pulse-energy-to-noise ratio, ^w^( t)dt/l!T^, can be of 
the order of 20 dB. 

By invoking the principle of superposition for optical j 
fibre transmission (as shown by Personick [5]) we can i 

extend our single pulse description to obtain a model for 
transmission of an entire pulse train. Thus, if we transmit^ 

j-: 

a sequence of on or off pulses, the photo-detector output ; 
r(t), can be written as 




r(t) = w(t-tj^) + n(t) 

k 


( 2 . 1 ) 
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where the time instants tj^ form a Poisson process having 
intensity "with 


^ (t) = 2 Oi P(t-IT) +% 


( 2 . 2 ) 


where 


P(t) ^ o 


is pulse response of fibre 


( 2 . 3 ) 


=0,1 are equiprobable data symbols 

\ Q ^ ^ dark current 

T is signalling interval 

n(t) is an additive noise process 

gjj. is avalanche gain factor 

w(t) is output pulse of photodetector. 

The number k of photon arrivals is random whose probability 
distribution is given by 


K 


prob(k = K) = 


A 


-A 


El 


( 2 . 4 ) 


where 



A (t )dt 


( 2 . 5 ) 


is the average number of photon arrival in the interval 
0 to T. . 


The subject of our investigation is : How should the 
photo-detector output r(t) = I(t) + n(t) be processed so as 
to minimize the probability of error? 



10 


Before answering this question in the next chapter, 
we briefly describe some of the previous work in the 
related field. 

2.2 PREVIOUS RELATED WORK 

In the past few years a useful communication theoretic 
model for data transmission over the fibre optical channels 
has evolved. Bar-DavM [6] and Gjvgliardi and Karp [7] 
ha,ve considered the optimal reception problem in the 
absence of dispersion (Inter symbol interference) and 
additive noise. Bar-David has considered the situation 
when the additive noise is also not present. The maximum 
likelihood function derived by him requires the knowledge 
of the electron emission instants. The receiver obtained 
is a weighted photon counter. If the output pulses are 
narrow enough to be considered as delta functions then 
this can be replaced by a linear filter with impulse 
response. 

h(-t) = In ( 

and optimum threshold is given as 

F = dt - dt (2.7) 

Gagliardi and Karp have considered the problem of 
M~ary Poisson detection in a back ground of additive 
Poisson noise. 




( 2 . 6 ) 
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Pei’sonick [8] has considered the detection problem 
in the presence of thermal noise and intersymbol inter~ 
ference. He presents a linear sub-optimxim solution. He 
gives the design of a front-end amplifier to take into 
account precisely the input pulse shape and equalizer 
filter shape. The output pulse shape after equalization 
is assumed to be a full raised -cosine waveform. 

Other sub-optimal receivers in the presence of noise 
and ISI are considered by Messerschmitt [9]- The receiver 
minimizes the mean square error (MSB), both for linear as 
well as decision feedback equalizers. The receiver is 
derived assuming the structure to be a linear filter 
followed by a forward and a backward tapped delay line. 
Receiver performance in terms of probability of error has 
not been evaluated < 

The MI detection has been considered by Foschini et al 
[l0]. The receiver obtained consists of two Viterbi- 
decoders in cascaiie-- Some of the relevant details of 
this receiver are presented in the next chapter where a 
modification of the same is considered. 



CHAPTER 3 


OPTIMJUj AMD SQB-OPTIfrlAL PROCESSING OE RECEIVED 

SIGNAL 

In this chapter we begin to answer the question posed 
at the end of section 2.1 (How to optimally process the 
phdtode tector output?) by presenting a derivation of the 
likelihood function associated with the received signal. The 
likelihood function so obtained is too complicated to yield 
a useful receiver. Under the assumptions of a high signal 
to noise ratio as obtains when an avalanche photodiode or 
a photo -multiplier tube is used» and quantized time 
instants of electron emissions, following {j^ » show that 
the optimal receiver consists of two Viterbi decoders in 
cascade. As hardware implementation of a Viterbi decoder at 
high data rates, with existing technologies^ is nearly ruled 
out, in this thesis we consider a modification of the likeli- 
hood function to obtain a decision feedback equalizer. 

5.1 GENERALIZED MXimM LIKELIHOOD DECODING 

In the digital system of Eig. 2.1 , the decoder has the 
task of observing the detector output and decoding the trans- 
mitted signal. The presence of added noise fields during 
photodetection, along with the added system noise prior to 
decoding circuitry, i causes errors in the bit decisions made. 
An accepted procedure for decoder design for generating these 
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decisions makes use of the concept of Maximum Likelihood (ML) 
test. Under this criterion, the decoder upon observing the 
noisy detectOT output, selects the bits that have the highest 
probability of causing this output. Mathematically, this 
requires the decoder to compute conditional probabilities 
p(r(t)ji) (called likelihood functions) for each combination 
of bits over TT sec ( 'C sec. i,s . observation interval), and 
select the maximum. The decoder computes 

= P(3:^(t)ji) = ike probabilitsf of receiving r(t) given 
^ ith intensity as received (3.1) 

for each i and selects the maximum. 

I 

For equally likely intensities the ML test is equivalent 
to MAP (maximum ■ aposteriori) test in which the detector 
selects the intensity^ which has the highest probability of 
having been sent^as the received intensity. A MAP detector is 
a minimum probability of error detector. Since for equally 
likely intensities, ML .teat is same as MAP test, it is also 
a minimum probability of error criterioh* 

If the received signal intensity happens to depend upon 
a random parameter 6, havifig apriori density p(0), then the 
■previous discussion is still valid, provided that we remember 
that p(r(t))i) is a conditional average over 0 [p] . That is 
we use p(i*(t)|i) = Bq p(3^(t)|i,©) » where p(r(t)|i,©) is 
the probability conditioned on the ith intensity being 
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received with parameter value ©« Thus likelihood function 
he c ome s 

\ ~ J^P(r)(t)| i.e) p(0) d© (3.2) 

This is called a generalized likelihood function. This 
concept can he generalized for several random parameters. 

3.2' GENERALIZED LIKELIHOOD EUNCTION EOR OPTICAL EIBRE 
DETECTOR 

If we denote the received signal hy 

r(t) = lj±) + n(t) (3.3) 

where is the information carrying filtered Poisson 

process 

k(t) 

IjjiCt) = -S' gj w(t-t .) (3.4) 

j=1 

and where index m (corresponding to \^{t) is hidden in the 
statistics of ^tj^l and k(t). These statistics are described 
hy equations (2.2) to (2.4) with A(t) replaced hy . A. m ^ ^ * 

As we have mentioned earlier^ the random variables 
represent the random gains g. w(t) is the pulse shape due to 
emission of an electron. In actual practice the additive 
noise at the output of the detector is not white hut it can 
he whitened hy a filter before additional processing. The 
effect of this filter and pulse shape due to emission, of an 
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electron can be combined together, then w(t) is this resultant 
pulse shape . 

The likelihood function conditioned on 
standard form [2] 

X 

I„(<:)r(t)dt - I^(t)dt) (3.5) 

0 0 o 

The desired likelihood function is the expectation of (3«5) 

with respect to I (t) for fixed m, i.e. 

m 

l^(r(t)) = ^ ( 3 . 6 ) 


Once the intensity ^(t) is specified, the above expectation 
is taken with respect to the number of arrivals, the arrival 
times and the avalanche gain values. The detailed evaluation 
of this expectation and the interpretation of the resulting 
structure in terms of simple physical operations on r(t), are 
our objectilres. The exact structure turns out to be so complex 
that many judicious approximations will have to be made to 
glean the essential nature of the operations. 


Neglecting the edge effects on the integrals and assuming 
that the observation time tT is much larger than the effective 
duration of a single pulse w(t), we can express the inner 
product and square term indicated in (3*5) as y 


[ Ijj^(t) r(t) dt = ^ gj P(t^) 


0 


(3.7) 
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where 


w(t-t^) r{t) d(.t 


(5.8) 


The square term is written as 
r k k 

= 2r .21 SzS-i E-( i ) 
0 “ 0 ^} 1=1 3 1 3 1 

where 


TT 

R(t) = j"w(u) w( t -u)du 


is the pulse correlati''a function. 


(5.9) 


(3.10) 


Substituting (3*7) and (3.9) into (3*5) and (3.6), we obtain 

k 

\(r(t)) = E^f exp(^- Z. ^^^3^ “■ 2^ 2 


L 




o 3=1 


o ^3 % 

2 fi-— 1 

R(t.-ti))^ (3.11) 


Employing the vector notation 

Sjj = (g^»g2» §15-) and 


~k “ (‘^■j»'^2’ "^k^ gives the expression 

k 

1 (r(t)) = Et e k rexp(|^Z g. R(t.) - 

Ir 2“^ ' 


N Z ^1 B.(t .~t^)) 


2N ^ “3 

° 3,1 


( 3 . 12 )‘ 


The joint density p(t^ ♦t 2 , ...» given as (see [f) , 

p. 1 1 1 ) 

n 

exp(- jA(t)dt) 


n' 


3=1 


(3.13) 


n 


After performing the expectation in (5*1 2) using (3»13)» 
we obtain for the likelihood function 


L (r(,t) = exp(-A (t)) 

“ n=0 o 


2:, y2Lj 


n 


j=1 
n 

exp(~~ 2 g. P(t ) - 21 g. gn ECt.-^t )) (3.14) 

^0 3=1 ^ ^ 3,1 ^ ^ 


n 


where p^^Cg^) is the discrete probability density function of 
the avalanche gain and where it is understood that when n = 0, 
the sumand is tasen to be unity . 


5y suitably normalizing the likelihood function, (3.14) » 

^ can be replaced by the pulse signal- to- noise ratio. The 
0 

normalization entails replacing E(t) by R(t)/R(0), P(t-) by 

J 

P(t^)/R(0) g and gj^ by g^^/g; consequently 

= g^ R(0)/R^ ( 3 . 15 ) 

and may be viewed as an average pulse signal-to-noise ratio. 

2 

Depending upon the presence or absence of avalanche gain, S 
is either large or small^ To save on notation we still repre- 


sent R(t)/R(0), gjj./® P(t^)/R(0)g as R*(t), g^ and R(t^) 

2 

respectively. Presence of S in expression will indicate that 
we are using the normalised quantities. Now the likelihood 
function becomes 


L^(r(t)) 


n 


= exp(-A(r)2 Ps(ej) , 

^ O n 

. .exp^S^( ^ gj ■■ 


(3.16) 
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We see that the likelihood fimcti on given in (3*16) 

is too complicated to be interpreted and to be implemented. 

In the presence of avalanche gains, the average signal to 
2 

noise ratio, S , is large. Ihis implies that the photon 
arrival times can be accurately estimated and these estimates 
can then be used to aid the detector in making accurate deci- 
sions* Heuristically , the receiver attempts to ' ’Vhiten" 
or peak up the pulse wit). The presence of G-aussian noise, 
however small, prevents pulse whitening via linear filtering. 
The nonlinear manner in which the receiver estimates the 
arrival times is of our interest and will be presented in 
the sequel. 

3.3 THE ASYMPTOTIC (8^— >00) LIKELIHOOD EUECTIOH 

In this section, the' basic idea is to asymptotically 

2 

evaluate the multiple sums or integrals. When S — >00 , the 
2n-fold integrals appearing in the likelihood function (3*16) 
become increasingly sensitive to the value of the exponent, 
and in the limit the integral is completely determined by 
the coordinates that maximize the exponent. This statement 
is made precise by the multi-dimensional version of Laplace 
Theorem [l lj which gives for each n 
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lim I dtn ,2, 




it «3 - 

n 

yz 21 6j S]^ R(tj-t]_) )] 



where (t*, ••., t*) and 


(g-, , g2* •••» maximif^e the exponent 

2«j p(tj) -yaz:®^ «iH(tj-tp 

D=1 j,l 


(3.17) 


(3.18) 


under the constraints o t^ ^ 71 , i = 1,2, . .. , n. 

and £g 

appear very difficult. For the time being let us assume 
that these sets are determined somehow. Rewriting right hand 
side of (5.17) as 

it Ps(s*)A(t*) exp[ s2(Z: g*J(tp - 

3=1 n 3=1 

'/2 2 Si - ‘^1^}]= exp(g^ B^(t)) (5.19) 

3»1 

where we have indicated the dependence of both the coefficients 
and exponent on the observation interval . Using (5*16), 
(3.17) and (3*19) we get the Dirichlet series 


The determination of extremiging values of 
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/ / \ \ jrs. O 

exp(-. 3 X (t)dt) y' G (t) exp(S^ B ('^) ) (3.20) 

n=0 

As S-^cO, it is -well known that the Dirichlet series is dominated 
by the term with the largest exponent, i*e., 

lim L (r(t)) /X-' exp(-A^) C*('C) exp( (3.21) 

3 CO ^ m n n 

"X" 

where B^('t-') is the largest exponent. ¥e see that n* is an 

estimate of the number of Poisson events occuring in the 

interval and that t*, t^j ...» t** are estimates of these 

occurrence times, while g*, g^, ...» g** are estimates of 

avalanche gains. Since neither the exponent in (3.21) nor 
n* 

the term 'J~\ pQ(gj) in Cjj^('tO is hypothesis sensitive, 

0=1 

the relevant portion of the likelihood function is 

n* 

\(r(t)) /-V/ exp(-A^) n (3.22) 

0=1 

where n* is the number of time points that maximize the 
exponent of (3.20) and .?t. are the values of these time 
points. We note here that likelihood function given in (3.22) 
is same as, the likelihood function derived by Bar-David j 
where he assumes that arrival times of photon are known exactly. 

Once the exponent is jointly optimized with respect to 
t^ and the estimates of the avalanche gain is not utilized 
further. This is so because the avalanche gain is a property 
of the photo-diode and conveys no information concerning the 
intensities. 
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Aa we have seen earlier that evaluation of (3*1 7) is 
a difficult task. It becomes simpler if we assume that the 
electron emissions are discrete. This we do in the next 
section. 

g.4 LIKELIHOOD EUDCTIOH WHEN ELECTRON EICLSSION TIMES ARE 
DISCRETE 


We assume that photons are constrained to arrive at 

discrete instants jA (1=1 »2, J = . This gives 

rise to the discrete likelihood function 

06 



where take value jA, (j = 1,2, ...» J), Using 

the same reasoning as before we will get expressions similar 
to (3.17) — (3.22), where expressions will contain sumnations 

instead of integrals. We choose the time quantization interval 
A, small enough to ensure that the probability of more than 
one electron emission occuring in a time interval A. is 
vanishingly small under each hypothesis Am(t) (effect of 
varying A is discussed in next chapter, Sec. 4.5). In this 
frame-work, the set of time points ^t^^J , specifies J-points 
in the interval (o,71 ), and the exponent can be written as 



2 Sj Ij - ’/2 S 

i=1 i,l=1 


g3_ <l3_ R(^3 - ^1) 

(3.24) 


where q.^ = 0,1 depending upon the photon arrival interval j A . 
It is clear that the product g. q. is inseperahle in the 

J J 

optimization of (3.24). Once the optimum values of g. and q. 

J J 

are determined, only the values of q. play a further role in 
the detection procedure. Keeping this in mind we let 
hj = gj q^t where will range over the allowable values of 
as well as zero. Por convenience sahe we call this discrete 
set B. Then, the optimization problem posed in (3.18) and 
(3.21), „ 

max (2. g P(t ) -'12 2 g, g, E(t.-t,)) (3.25) 

K .t^,n 3=:1 3 3 0 1 3 J. 


becomes 

J J 

max ^ b. P(jA) -y2 . b . b, B.( 3 ,<h -Iji ) 

b, ,b« ... b, • ^ ^ 3»1=1 ^ (3.26) 

J 3=1 


Here we note that maximization in (3.25) requires finding of 

i 

n J 

maximum amongst these maxima. Such a complicated procedure 
is obviated in (3.26) by fixing J. We recognize the expre- 
ssion in (3.26) as the familiar likelihood function for 



and 


i_q 


for each n and then choosing the 
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estimating in presence of Gaussian noise when these are 
transmitted through a channel with impulse response w(t)» 

Thus can be estimated using Viterbi algorithm as des- 
cribed by Forney [l 2j . Now the likelihood function can be 
written as 

\(r(t))^'^ expC-Ajjj) ^ ^ (3.27) 

Substituting X(t) from (2.2) in (3.27) we get 
li (r(t)) ^ exp(- / (2 a^ p(t-nT)+>^Q)dt) . 

* TT ^ (5.28) 

3=1 

N 

Now the data sequence ^a 1 can be estimated by calculating 

u a=1 N 

the likelihood function (3.28) for all the possible 2 sequences. : 
Since p(t) in (2.2) has finite duration, the likelihood function 
in (3.28) can be calculated recursively and can be 

maximized using the Viterbi algorithm. Then the receiver stru- 
cture is as shown in Fig. 3.1 . 

As we know that implementation of Viterbi-algorithm is , 

quite complicated, in the next section, we derive a sub-optimum i 

structure by assuming that past symbols are correctly known to ; 
the receiver. This gives rise to the decision feedback equalized. ' 



Fig. 3.1 FIBRE OPTIC RECEIVER 
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3.5 DECISION FEEDBACK EQUALIZER 


Ihe optimization of equation (3-26) can be done by 
Viterbi algorithm as pointed out earlier. ¥e also know that 
a sub--optimum solution of this problem using decision feedback 
equalization exists. This is given by Austin |j 3 J and Salz 
;£4j . Hence we do not dwell upon optimization of (3.26) and 
assume that above receiver structures can be easily implemented. 
Now we go on to optimization of (3*28), Let us denote the 
impulse response p(t) of the fibre as 


(t) = 


/p(t) 0 :< t < I 


otherwise 


(5.29) 


P2(t) 


p(t+T) 0 ^t<I 
otherwise 


xP(t 

L o 


p (t) 


p(t+{m~1 )T) 
0 


0 < t < I 

otherwise 


where 

p(t) is non zero over the interval 0 to ml and zero for t 
When we are making decisions bit by bit then it is easy to 
note that if we are estimating the kth bitj, our observation 
interval in (3.28) should be (k-l)I to (k+m-1)I. Keeping this 
in niiiid W6 can write (3*28) as 
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(k+m-1 )T 


L (r(t)) = exp(- ^ ^ p(t“nT)dt+A 

“-'(k-Di 


kj k-1 

•IT ( Z: 

J=(k-1 )J+1 “=k-®+1 “ 


(DA-(k-l)jA) + 


\V^i2A- (k-l)J^) +k^) 3 
(k+1)J k-1 

• TT (Z. a„s( 0/Si-kJA)+a. P-C DA~kJA) + 

3=kJ+1 • n=k-mH-2 ^ k-n+2 ^ ^ ^ ^ 

kJA) +Aq) 3 


(k+m-1 )J 

I I ( 3^Pjjj( “(k-Hn— 2) J^) + 

j-(k-Hn-2)+1 

k-fm- 1 q . 

2 Vk+m-a‘ +\) ’ <5-50) 

n=k+1 

Here we have split the terms in interval (k-1)T to kT, 
kl to (k+1)T, •..» (k4m-2)T to (k+m-l)!* At this point 
the natural thing to do would he to take the expectation 
of the likelihood f-unction in (3»30), i.e,, ]j(r(t) a. ) = 

K f 3,-“ 

(L(r(t)|a. -. +))• The likelihood function has the 

form; Zj exp(- /X (t) dt Jf rXjj.( D^) Here n take 

2^ values. Implementation of this likelihood function is 

difficult because we have to perform multiplications in 2^ 

channels and then add the results. Also in each channel the 

number of multiplications keep changing from interval to 

interval as the q.’s which are non-zero keep changing in 

J 


s 
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each interval. This likelihood function is to be calculated 
for each hypothesis (a^ = 0, a^j. = 1)* take the 

expectation of log- likelihood function -with respect to 
future symbols, the receiver structure is greatly 
simplified. This gives a structure which is not optimum, 
but it is simpler to impletoent. Taking log of (3-30) 
them averaging over eqUiprobable future symbols we get 

2! % \+1-n " ®k^1 

* ^ n=k^m+1 

j ^ 


k-1 


’ 2^ ®'n ^k+2-n ' ®k^2 ” 


G=k-m+2 ^ 

+ y2 Vk+2-n^^^'^'^)^ 

^k 


ln( 2] \+2-n ( 3^) +p^ ( 3 A) ) 

k-1 

21 S\+r-n " \V^2(P^+P2+...+Py„^)-Xo' 


n«k-m+r 

+ ^ ^ «3+(k4.r-2)j[ ^ 

2 k 
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©1 »©2****^°r-1 
c^=(0,l) 


kktZ „ 5k+r-n'3^>+'\= + 

n;=k--m+r 


2 ©i Pi( 3^))J 


contd. . . 
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••• +P„.,) 


1 


j=1 


m-i c^ 'i3+(k-Ha-2)Ji 


Ao+ kl. c^P^(3A)J 




r"2^ 

C^=(0,1) 


jr C i 

m-1 




(3.31) 


where we have used 


r 

j K(t)it = P 


.T 

I 

'0 '3 3 

A zero or one is decided depending upon 


(3.32) 


(3.33) 


Substituting (3.31) in (3.33) and cancelling out the common 
terms we get 


m J 

"5^ -D 5*' 


+ y2 2 ln(1+ezj 


ln(1 + 


k-1 


p^(da) 




) 


h=K:-m+1 


P2( 3A ) 


) . 


Z: Vk+2-n‘J-)-^^o 

nek-m+2 


(1 + 


S Vk+2-n'3^> + + 

n=k«m+2 


P2(3A) 


) 
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J 


j=1 


^3+(k+r-2)J 


z: 

c^=(0,1) 


ln(l + ^ 


Pj,(3^) 


.) 


n^-m+r ^ 


-£ZT ^'^'lj+(k-Hn~2)J ZL ln{l+ ^ 


+ 


c^=(0,1) 


(3.34) 


i 

The receiver so obtained is shown in'Pig- 3.2. In the next 
section we discuss some of the hardware features and com- 
plexity of the receiver. 

5.6 RECEIVER OOMPI^JXlTTr 

The receiver shown in Fig. 5-^2, requires a memory of 
size 2'^’*’lcJxot where m is the number of bauds for which impulse 
response r(t) is non-zero, J is no. of samples per baud 
CJ = T/^ ) . Since we are required to store values 
ln(1+p^(3/6v)) in the memory^ we give some estimate of word 
size. For the impulse response considered (shown in Fig. 2.2) 
the maximum value is ln( 1+37.3) ■•=^.7 and the minimum value 
is ln(1+2.1x10'’^) = 2.1x10“^. For this range of numbers 
we require 12-'bit words. The (m-1) serial shift registers 




JQDmBDVFlBHB 
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required to store q.'s should .be J bits long. each. Then we 

3 

require m-summers and m-acounrulator registers. Since 
average number of photons received is 100 to 200^ the threshold 
will be around these numbers. Hence accumulator register 
should be 17-bits long. ¥e will show in the next chapter 
the number of samples required per baud is approximately 
twice the number of average electrons » hcncO| J will range 
from 2.00 to 400. 



OHAPCER 4 


PERPOBj^IMGE OR SHE DECISION EEEDBAC3K EQUALIZER 

In, this chapter we evaluate the performance of the 
Decision 1'eed‘back Equalizer obtained in the last chapter* 

Also, we compare its performance with an ’Integrate and 
Dump* receiver, ’One Band’ receiver and ‘Ideal’ receiver* 
Numerical calculation of the probability of error calls for 
a great deal of computer time (approx. 30 min.), over- 
come this difficulty, we calculate an upper bound on the 
probability of error. The upper bound on the probability 
of error is obtained by considering the situation in which 
intersymbol interference effects are the worst that can 
occur. 

The inter symbol interference may aid or hinder the 
detection of a symbol depending upon whether the symbol is 
'1' or ’O’. When we are trying to detect a symbol ’1* in 
the prcvsence of ISI, and the symbols preceding and follow- 
ing it are all M', then these will help in detection by 
producing more electrons. But if we are detecting a symbol 
’O' in the above situation, the probability of error will 
be large. Similarly, detection of a ’1’ in the midst of 
all zero's will have a probability of error larger than 
that for any other combination of the preceding and 

following symbols. Hence two worst conditions 

CENTRAL UBRAIlf 
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detection of a ’1’ in all zero’s, and dcteotion of a zero 
in all one’s, These give the upper bound on the probability 
of error. In calculating these probabilities, 'we fixed 
the threshold such that the probability of error was same in 
both of these worst cases. 

Closed form expressions for the probability of error in 
the above cases are not available. The probability of error 
can be evaluated only by numerical methods. One such method 
suitable for our purpose is Gram-Gharlier series expansion 
method We will describe this method In the next 

section. 


We have also compared our analytical results with 
simulation studies. Since large computer time (approx 30 rain) 
is required, we simulated the receiver only for the case when 
Aa 50 (« average number of electrons emitted by a single 
pulse). 


4,1 GEAM-CHAELm SERIES EXPANSION METHOD 



h(t) 


x(t) 


Elg. 4.1 




Let be a series of impulses whose arrival times are 
Poisson distributed with the mean intensity The moment 
generating function ©f the randoia variable x(t) is 
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M^(s) « E(exp(sx('t)) 


k(t ) 

E(exp(s h(t-t^))) 


(4.1) 


where k(t) is the number of counts in the interval (o,t) and 
are arrival instants, To calculate this expectation, 
first we find the conditional expectation conditioned on the 
number of counts, and then average over the number of counts. 


i » e • , 


:^(8) s E^PE^(exp(s 2" J 

^ *** J h(t-t^) p(t-j jtg,. .. ,t^|k) 

^0 0 o 

at^ atg ... dtj 


M. 


(4.2) 


B,, 


k t 


5 exp(s h(t-t.) -- r ^- dt J 


where 


* 2 

t 

:^(s) s ^ exp(s h(t-.z)) 


dz 


A 


(4.3) 


Averaging over k yields : 

K 


M. 


oo A 




fcaO 


(4.4) 


exp Rm ( s ) '• 1 ) 

cf (exp(8h(t-2))'>(z)->(2))4* j 


«» exp 
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Expanding thd inner exponential term fey Taylor series » we 
get j 

Oo t 

M (b) w exp(2r fr C fe^(t««) X(z)dz) (4*5) 

^ Q=1 

0 

Therefore, this has cummulants (semi-invariants) as 

y'q ~ ^ lX(z)dz (4.6) 

0 


Since the numerical values of these cummulants' increase with 
q, we consider the normalized process 



as 


x(t) - yc^ 



(4.7) 


This has moment generating function as 


M^(fl) « E [^exp(sx)3 




CO 


Z oT yTqZa ’q ■ 

*1 2 


q= 

CO 


iZl' 


exp (2 

q^l 


ql 


(4.8) 


(4.9) 

I 
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■where 


r, 

r^ 

9. 




Q >2 


(4.9) 


One ■way to obtain the density function of a random 
variable from knowledge of its commulants is to use the 
G-ram-Charlier series expansion of the density function, 
namely 


oo 


p(x) =2- ^n 
n^O ^ 


(4.10) 


where 


The corresponding moment generating function of p(») is 

60 

M^(^) * J exp(sx) p(x) dx 

OO ^ n 

« 51 C ! exp(sx) -==: exp(-x /2) dx 
“ n J dx*" 


^<p 


n«sO 

y s“ exp(s^/2) 

^0 

£/i ^ «P(»V2) 

n»0 


(4.12) 


where 


/*n - °n 


nl 


(4.13) 
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The coefficients, G^, in the expansion (4.10) can be 
calculated from knowledge of oumraulants, 'by equating the 
moment generating function ealculated in (4*8) and (4*12), 

2 jT" -h o w A 

expCf- + 2,^ fr) * €xp(s^/2) 2 )U,14) 

n=3 * nwo 

Expanding the left-hand side expression, after jjpnfelling 

2 

the common factor exp(s /2), we get 



Equating coefficients of both sides of (4*14) using (4,15) we 
get 

yu . 1 (4,16) 

A^1 “ ® 

^2 =» 0 

m -2 

The coefficients, G^, can be calculated using equation (4. 13), 
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3!he functions* can alsc la« calculated recursi- 

vely by noting that 

H^(x) « 

where 

Hq(x) *» 

and 

H^(x) w - ~5~expC-x^/2) 




expC'-x / 2) , 


The probability of error can be talculated uains 
followii^g two relations 


X. 


r /• <» 

3 p(x)dx ^ ^ C^H^(x)dx 

•oo 


»-«c n«o 

X. 


CO 


X. 


S H,(x)d* + 21 ®n S„ ®^'*>** 

«> -cc 

<X? X^ 

5 4r •xp(-x*/2)ta ♦ X 

« 4^^ b-5 '* 


00 


4 erro(-Xg) *o<' ® 

a»5 

(4.18) 


y 

where erlb (y) *= ‘^’•C 


5 J 7X ®3tp(-3c /2)dx 
0 


(4.19) 


Similarly, we have 
00 

vf p(x)dx « i erfcCx )- (4.20) 
x^ n»'5 


00 
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The equalizer v« derived in the last chapter (fig*3*2) 
operates the estimates of arrival times* But some 

of these estimates will he inoorrect because of the thermal 
noise in the first detector. So we can write the. output of 
first detector as 

= q + p . (4, 21 ) 

w j 3 

where * p. =: 1 means we detect an electron 

J tJ 

emission when in fact there is none* 

|( Similarly, p^ » -1 means we fail to detect an ©aitted electron, 

p. = 0 means no error. To take into account this situation, 

J 

we assume that input to filter^ h(t), in 3'ig. 4.1 is 
lc(t) k^(t) 

2 ^ ,^(ts.t^) + ^ (4.22) 

3=0 

The first terra is the same as in the previous discussion* 

The second tem takes into account the errors, 

is arrival time of those error impulses, and k^(t) is 
the number of errors in the interval o to t. 

We now make some simplifying assumptions about the statis- 
tics of the error process. We assume that 

(1) are independent of each other and 1 are 

equiprobable. 
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(2) Ihe arrival time 2^ is uniform In time interval o-T. 

(3) The error process is independent of the process 

The assumptions (1) and (2) are not quite correct 
because the errors in the first DSTl -will aetuallj be in 
bursts. But for simplicity we assume^the /errors to be 
random. Then, by assumption (1), the number of errors in 
time T will obey binomial distribution, which can be 
approximated by a Boisson distribution. Hence we make the 
fourth assumptions that 

(4) The number of errors, k^CT), in time interval of T 
is Poisson with mean e« 

Now the output of the filter h(t) due to the error 
process ^ 2^^ will be 

kg(T) 

B(T) « ^ h(T-2^) ^^*23) 

jssO 

Proceeding as before, we can write the moment generating 
function of the output error process as 


I^(e) , B [^xp(s B(T) ^ 

■ ■«. I exp(se^h(T-Z) ^ J 


(4*24) 


j=1 0 


Since e, » i 1 equiprobable we get 
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Mj.(s) 

= E, 
k 

\ i ($) 

2 r«o 

■where 

A 

T 

- i 

j exp(sh(!r«.z).)az 


0 


and 

T 

^ exp(-s h(!I!-z))dz ■ 

B 



0 


Simplifying 

(4.25) we get 

M5,(s) 

= \ 

/ A-f-B vk*^ 

I 


CO 

% 1 

k=o 

|r exp(-e) (^)^ 


= exp 

[^(A+B . .,,3] 



T 

PV* 


= exp 

1 r exp(sh(T-z)dz 



0 


(4.25) 


(4.26) 


(4.27) 


T 


Expanding the inner exponential term, v/e get [. 

oO CO T I 

M^(s) « expPfX ^ C h'^d-z) f dz 4 f ^ i?(t^z)dzl 

^ q=1 qi q=1 Q. 


q=1 q 

vy 

OO 

“ 2 - IT /'q 

q=2,even 


A 


(4.28) 
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■where ^ 

; q even (4.29) 

0 

Since the processes and ^ are indepoadent, the 

moment generating function of the output process will he obtained 
by multiplying (4.28) and (4.5). 


M^(s) = M^(s).Me(s) 



q 

Z ^ /3) 

q=2,even 


(4.30) 


The expressions (4.7) to (4.16) are still valid if we replace 
T^q where it is understood that *= 0 for q odd. 

Hence, we can calculate the probability of error using (4.18) 
and (4.20). 

4.2 HUMERI Oiii RESTJMS 

In the previous section we outlined a procedure to calculate 
the probability of error. In this section we use this procedure 
to calculate an upper-bound on the probability of error by 
calculating the PE for worst case of ISI. ¥e compare the 
performance of the decision feedback equalizer with several other 
receiver structures. These receiver structures are 

(1) Integrator and Pump t Here the filter h(t) has the response 
h(t) w 1; o<t<T, and h(t) »» o otherwise. This is the simplest 
receiver and in fact, is nothing more than a counter. 
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(2) One BsWd receiver ; This is also a decision feedback 
equalizer. But here "we decide about the symbol after the 
first baud only. Therefore, there is no future symbol inter- 
ference# But there is a loss of signal energy as the energy 
in the first baud only is being utilized. 

(5) Ideal : In this receiver we assume that only the symbol 
being detected is unknown to the receiver, all other symbols, 
past or future are known to the receiver. 

Por the purpose of calculations we have used the impulse 

response of fibre as shown in Pig. 2.2 with Z » 10 km. The 

time interval T is taken to be 150 ns. When an average of 50 

primaiy electrons are emitted by the photo -detect or, then the 

number of electrons in the first and second baud are 27.7 and 

20.1, respectively. We assume that the response extends upto 

three bauds only. The upper baund is calculated by calculating 

the error probability for the sequence 00100. The coefficients, 

0 of the G-ram-Oharlier series were obtained by eqns. (4.6), 
n 

(4.9), (4.15) and (4.16). Per calculating the cummulants 7q 
the appropriate impulse responses h(t) were taken. These impulse 
responses are i 

Per ’Integrate and Dump’ receiver 

h(T-t) * 1 0 ^ 't < T 


0 


otherwise. 


( 4 . 51 ) 
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Por ’One Baud’ receiver 


(4.32) 

Bor ’Decision Beedback Equalizer’ 

Pi(t) 

h((k+2)[ll-t) = ln(1+ 7r-rrri -tt t -t ' " 

P2(t) P2('fc) 


4;(lc-1)Ti$t<kT 


_. i 

— 4 , 


kT 4t < (k+1 )B 

, V-rit) 

1„(U _)(u^_f^)(u 


PjCt) 




Bor ’Ideal’ Receiver 


Pi (■fe) 

h((ic.2)i-t) = in(u M 

Po(t) 

}; (k!E)^t <(k+1 )T 


= In ( 1 ■<• 


p,(t) 


Ihe probability of error -was calculated using cqns. (4.13) to 
(4.20). 

lable 4,1 shows the probability of error for different 
types of receiver when there Is no error in the estimates of 
arrival time of. photons, The probability of error is plotted 
in Big. 4,2, In fables 4*2 and 4,3 the probability of error is 
calculated when the average number of errors in the estimates of 



Table 4.1 

Probability of error when 
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arrival time is 2?^, and 5^, of the average number of photons 
received* IThe results are plotted in Pigs. 4.3 and 4,4. 

Prom the plots in Pig. 4.2, 4.3 and 4.4, we can conclude 
that DPE is 3.5- dB better than as Simple Integrate and Dump 
receiver. It is better than a One Baud receiver by 0,S2 dB. 
But the DPE is inferior to Ideal receive by 1.4 dB. These 
calculations are for the case when the probability of error is 



To see the effect of errors in the estimation of arrival 

! 

time, we have calculated the probability of error for different 
values of e/^ (e * average number of errors in one baud). The 
results are tabulated in Table 4.4 and shown as a graph in Pig. 

4.5. Prom these we see that the system degradation is not 
very large, 

4.3 SIMULATION RESULTS 

The DPE in Pig. 5 . 2 was simulated for average number of primary! 
electrons equal to 50. Pirst, we simulated the receiver when the [ 
input sequence was 00100. The feedback filter contained 00 and 

the feedback path was broken. This corresponds to the model used ■ 

i 

in the previous section to calculate the upper bound oq PE. A t 
computer-run with 10^ iterations gave 2 -errors. Another run with li 
10^ iterations gave 18 errors. This gives the PE as 2x10 , 5 

times more than the value calculated in the previous section. The | 
time taken for 10^ iterations was approximately 30 min so we did , 


not make further trials 



Probability of error ’when (e//\) ='02 
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The simulation of actual BPS, when input symbols were also 

c 

random, gave 21 errors in 10 Iterations. This number is approxi 
mately same as above, indicating that there was no error propaga- 
tion in the PPE. 

¥e also simulated the One Baud receiver. Computer-runs of 

10^ iterations gave number of errors at 9,8 and 11 in three 

-3 

different trials. This is equivalent to PE = 10 , which is 

2.5 times the PB calculated numerically in the last section. The 
reason for discrepancy may be that v;e took the arrival times in 
the previous section to be continuous while here we are quanti- 
zing the arrival times. 

4.4 TIMING ERROR EPPEGTS 

In calculating the probability of error, we have been assu- 
ming that the bit timing is perfect and the decoder is perfectly 
synchronized. If timing offsets occur during a bit period, • 
because of inaccurate time synchronization of the transmitter and 
receiver, processing takes place in offset intervals, leading to 
system degradation. To get an idea of this, we assume that the, 
sequence OOlOO is repeated. Then we calculate the probability of 
error for various AT/T timing offset values. The results are 
tabulated in Table 4.5 and plotted in Pig. 4.6. The results 
show a large degradation (increase in probability of error) as 
the timing offset reaches *5. The timing offsets effectively 
decrease the signal energy end increase the noise. 
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Table 4.5 

Effect of time shift 



A= 100 
?e 

200 

Pe 

0 

0.2024430E-06 

0.2008073E-11 

.1 

0.1435179B-05 

0.1642984B-10 

.2 

0.1643697E-04 

0.1997740E-08 

.3 

0.16S9365E-03 

0.1 8911 OOE-06 

.4 

0.101 91 90E-02 

0.6147204E-05 

.5 

0.2832715E~02 

0.4366327E-04 

.6 

0.3594S22E-02 

0.61 8791 2E-04, 

.7 

0.1700239S-02 

0.1643831B-04 

.8 

0.3151896E-03 

0, 63285 24E-06 

,9 

0.1 6843 92E~04 

0. 20895 S6E-08 



ERROR WiT 
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4,5 EPPECT OP sampling RAPE 


The number of electrons emitted, n, by the photodetector 
is a random number, Poisson distributed with mean /s (average 
number of electrons, due to one pulse). Since we are taking only 
I(= P/^ ) samples per baud, we will not be able to detect all 
the emitted electrons. Phe probability that the number of emitted 
electrons, n is greater than J, is given by 

^ / s k’ 

P(„>j) . 21 C4.34) 

k=J+1 


Phis probability can be evaluated in terms of incomplete gramma 
function 3 

X 

^ exp(-t)t^”^ dt (4*35) 

o 

by using the identity 

k=J+1 


f 

i, 

t. 

[ 

E 

F 

I 


Phese values are tabulated in Table 4.6 for A= 50 and A® 100, 
Also, we have tabulated the average number of electrons that can 
be detected in one baud, Phese values are obtained by generating 
a Poisson random variable in the digital computer and then count- 
ing the electrons detected in a baud interval. The results are 
plotted in Pig. 4.7, "We see that if J = 2/\, the average, number 
of electrons detected is approximately A. i also the probability 
that n is greater than 2A is very small. So a sampling 



Table 4.6 

Variation in average number of Electrons 
emitted with sampling rate J 
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Sampling 

Average Fo. 

Prob(n? J) 

Average No. 

Prob (n > J) 

rate 

of detected 

/v 

of detected 

A = 100 

J 

electron n 

/\ 

electron n 



A = 50 


A = 100 


50 

36.39 

.546 



60 

39.44 

.104 



70 

41.37 

.534E-2 



80 

43.88 

.762E-4 



90 

46.98 

.335E-6 



100 

48.14 

.509B-9 

72.38 

.145 

110 

48.74 

,291E-12 

75.59 

.107 

120 

50.54 

. 686E-16 

79.38 

,120E-2 

130 

51.03 

. 7 IOE -20 

81.90 

.345E-4 

140 

50.95 

.333E-24 

34.13 

.478B-6 

150 ■ 

51.23 

.81 SB- 2 9 

86.37 

.333E-8 

160 



88.22 

.121E-10 

170 



90.66 

.242E-13 

180 



93.04 

,270E-16 

190 



94.16 

.174E-19 

200 


' 

96.74 

.619E-23 
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frequency •which is twice the number of average electrons 
emitted is sufficient. Moreover, it will not be required to 
detect the average of /\ electrons in every interval as this 
occurs only when the bit sequence is all 1’s. In all other bit 
combinations, the average number of electrons required to be 
detected will be less thanyA, 

The effect of the first DFE is to filter out the thermal 
noise. If J is large the number of errors e (errors in arrival 
time) will be large because in the first DPE the member of inter~ 
fering symbols will increase. Even if PE for this DIE remains 
constant, increase in J will mean a corresponding increase in 
e. An arbitrary increase in J defeats the purpose of the first 
EPS which effectively filters out the thermal noise. Hence, 

J = 2 A seems to be a reasonable choice. How actually the 
overall system performance changes with J can be answered 
only through simulation. Since simulation requires a large CPU 
time we have not done this but, have presented a beuristio 
argument , 



CHAPTER 5 


COHCIUSIOH 

In this thesis we have considered the problem of optimal 
processing of received signal in an optical digital fibre 
commiinication system. We evaluated the likelihood function of 
the received signal arid saw that the likelihood function is too 
complicated to yield a useful receiver. Thus, we made the 
assumptions that pulse-signal-to-noise ratio is very high (which 
holds when an avalanche diode is used) and electron emission 
times are quantized. The maximum likelihood function can then 
be mechanized using two "Viterbi decoders in cascade. The 
first Viterbi decoder optimizes the familiar problem of 
detection of symbols in Gaussian noise and ISI. Implementation 
of Viterbi decoders with the present technology is almost ruled 
out. So we replace the two cascaded Viterbi decoders by two 
decision feedback equalizers. The first decision feedback 
equalizer. is same as the receiver derived by Austin [13] and 
Salz C 14 ]. In deriving the second decision feedback equalizer, 
we averaged the log-likelihood function instead of the likelihood 
function. To evaluate the averaged likelihood function, multi- 
plications abe to be done in 2® parallel channels whereas only 
additions in ® channels are required to evaluate the averaged 
log-likelihdod function. Also, the memory size required reduces 
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from to 2’^'“'^xmxJ* Hence, we see that there is a 

considerable advantage in using the averaged log-likelihood 
function compared to the averaged likelihood function. 

In Chapter 4 we have evaluated the performance of the 
proposed decision feedback equalizer. Since exact calculation 
of probability of error is difficult, we have calculated an 
upper bound on the probability of error, using G-ram-Charlier 
series expansion -method. We saw that the performance of the 
decision feedback equalizer considered is superior to a simple 
'Integrate and Dump' receiver. Its performance is also better 
than a 'One Baud' receiver, but by less than 1 dB, Thus, when 
the impulse response of fibre is skewed to the first baud (as 
in the case of a short length fibre), the simpler 'One Baud’ 
receiver can be used without too much loss of performance, The 
'One Baud' receiver is actually a decision feedback equalizer 
in which ISI due to future symbols has been neglected. 

The receiver considered here has two main drawbacks. One, 
the sampling rate, J, required is very high, and the other, 
high sensitivity of receiver to timing offsets. The sampling 
rate required is twice the average number of electrons emitted. 
This high sampling rate puts an upper limit on data rate which 
can be transmitted for a given receiver complexity. Also, the 
memory required increases exponentially with the nUB|ber of inter- 
fering symbols. Hence the receiver derived is useful only when 



62 


the number of interfering symbols is not veiy large, 

There is a large degradation in receiver performance with 
timing offset errors. Hence good synchronization schemes are 
required for reliable communication. 

Some of the aspects of the present study which could be 
investigated further are techniques of optimizing equation 
(3,2-5) so that we do not have to make the assumption of 
quantized electron emission times. The second HHB will also 
change accordingly. The receiver proposed in this thesis has 
a drawback that sampling required is high which can be overcome 
by decreasing the sampling rate J stnd allowing q.j’s to take 
multiple values, instead of only binary values. This will 
result in complicating the first DFB and some loss of performence 
in the second DPE, 
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